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Abstract. We study stable solutions of the equation A 2 u = e u , both in entire 
space and in bounded domains. 

1. Introduction 

A classical result attributed^ to G.I. Barenblatt asserts that there exists infinitely 
many solutions to the equation 

(1.1) - Au = 2e u in ft, 

whenever ft is the unit ball of R 3 and the equation is supplemented with a homoge- 
neous Dirichlet boundary condition. The result appeared in a volume edited by I.M. 
Gel'fand [18], whose name the problem now bears. We refer the interested reader to 
the book [12) for some of the developments of this equation in the more than sixty 
years that separate us from Barenblatt's discovery. Let us simply mention that K. 
Nagasaki and T. Suzuki J^Hj complctly classified the solutions found by Barenblatt 
according to their Morse index. Much of what can be said of the equation posed in 
a general domain ft rests, through a blow-up analysis, upon Liouville-type theorems 
for finite Morse index solutions of the equation posed in entire space. This lead N. 
Dancer and A. Farina [9] to proving that in dimension 3 < N < 9, any solution to 
(jl.lj) in ft = M. N is unstable outside every compact set and so, it has infinite Morse 
index. 

In the present work, we consider the fourth-order analogue of the Gel'fand problem. 
Motivated by the aforementioned results and by the expanding literature on fourth- 
order equations, see in particular the books [17] by F. Gazzola, H.-Ch. Grunau, 
and G. Sweers, and [15] by P. Esposito, N. Ghoussoub, and Y. Guo, we want to 
classify solutions of 

(1.2) A 2 u = e u inl", 

which are stable (resp. stable outside a compact set), that is, solutions such that 

(1.3) / e u (p 2 dx< / \Aip\ 2 dx forallv3G J ff 2 (ft)n J ff n 1 (ft), 
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where ft is M^resp. the complement of some compact subset of R ). Consider 
first radial solutions. Noting that the equation is invariant under the scaling trans- 
formation 

(1.4) u\(x) = u{Xx) +41nA, xE~R N ,\>0, 
we may always assume that w(0) = 0. 

Proposition 1.1. Let 5 < N < 12. Assume u is a radial solution of (|1.2|) . Let 

v = —Am, = u(0), and f3 = v(0). There exist Pi > [3q > depending on N only 
such that 

• j8 > A>- 

• If f3 = /3q, u is unstable outside every compact set. 

• If /3 £ (/3o,/?i), u is unstable, but u is stable outside a compact set. 

• If ft > Pi, u is stable. 

Remark 1.2. The fact that no radial solution exists for (3 < (3q is due to G. Arioli, F. 
Gazzola and H.-C. Grunau [T|. In addition, E. Berchio, A. Farina, A. Ferrero, and 
F. Gazzola first proved in [3] that for 5 < N < 12, u is stable outside a compact seiQ 
if and only if j3 > (3q. The novelty here is the number f3±: our result characterizes 
stable radial solutions when 5 < A < 12. See [3], [32] for the remaining cases 
1 < N < 4 and N > 13. 

In particular, there is no hope of proving a result similar to that of Dancer and 
Farina in our context, without further restrictions on the solution. One might ask 
whether all stable solutions are radial, at least in dimension 5 < N < 12. This is 
still not the case. 

Theorem 1.3. Assume N > 5. Take a point x° = (x^, . . . , x%) € R^, parameters 
ai, . . . , apf > 1 + A/2, and let 

N 

(1.5) p{x)=Y. a ^- x ^ 2 - 

i=i 

Then, there exists a solution u of U.2\) such that 

(1.6) u(x) = -p{x) + 0(\x\ 4 - N ) as \x\^oo, 

In particular, u is stable outside a compact set (resp. stable, if min^i....^ &i is 
large enough) and u is not radial about any point if the coefficients ai are not all 
equal. 

Remark 1.4. Using the scaling (|1.4[) . one immediately obtains a solution u of (|1.2p 
such that u(x) = —p(x) + C + 0(\x\ 4 ~ N ) as |a;| — > oo, under the sole assumption 
that Oj > for alH = 1, . . . , A. 

Remark 1.5. In fact, any solution satisfying (|1.6[) has finite Morse index, thanks to 
the Cwikel-Lieb-Rozcnbljum formula (see G.V. Rozenbljum [26] and D. Levin [22] 
for the formula, as well as A. Farina |16] for its application to semilinear elliptic 
equations). We thank A. Farina for this observation. 



the notion of stability outside a compact set that we use here is stronger than the one given 
in [3]. One can easily check that the results of [3] remain true in our setting. 
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Remark 1.6. In dimension TV > 13, the radial solution upo(x) = —4 In \x\ + 0(1) is 
stable, see [3]. Observe that up does not satisfy f| 1 .6[) . In dimension 5 < N < 12, it 
would be interesting to determine whether, up to rescaling and rotation, all stable 
solutions do satisfy (|1.6[) . 

All stable solutions that we have encountered so far have quadratic behavior at 
infinity. In particular, letting 

v = —Au and v(r) = -r v da, 

JdB r 

these solutions satisfy v(oo) > 0, where 

t;(oo) := lim v(r). 

r— >+oo 

This motivates the following Liouville-type result. 

Theorem 1.7. Assume 5 < N < 12. Let u be a solution of (| 1 .2[) such that 
v(oo) = 0. Then, u is unstable outside every compact set. 

Remark 1.8. As observed in [3], in dimension N — 4, applying inequality (|1.3[) with 
a standard cut-off function, we easily see that if u is stable outside a compact set, 
then e" € L 1 (R 4 ). Thanks to the work of C.-S. Lin [23], up to a rotation of space, 
u must satisfy 

(1.7) u(x) = -p(x) - 47 In \x\ + c + £>(|.x|~ T ) as \x\ ->• 00, 

where p(x) is of the form (|1.5|) with a; > 0, 7 = J R 4 e w cte < 2, cq and r > 
are constants. Conversely, there exist solutions of the form (|1.7[) . as proved by J. 
Wei and D. Ye in [35] ■ If m addition zT(oo) = 0, then 7 = 2 and up to translation 
and the scaling (|1.4[) . 

u(x) = -4 In (1 + \x\ 2 ) + In 384, for all x g M 4 . 

Now, let us turn to bounded domains. We begin by recalling a few known results 
on the Gel'fand problem for the biharmonic operator, when the domain is the unit 
ball and the equation is supplemented with a homogeneous Dirichlet boundary 
condition, i.e. for A > 0, we consider the equation 

A 2 u = Ae" in B, 

u = |Vit| = on dB. 

It is known that there exists an extremal parameter A* = A* (TV) > such that 
the problem has at least one solution (which is stable) for A < A*, a weak stable 
solution u* for A = A* and no solution for A > A* (sec G. Arioli, F. Gazzola, F.-C. 
Grunau and E. Mitidieri [2]). The unique weak stable solution u* associated to 
A = A* is classical if and only if 1 < N < 12 (see the work by J. Davila, I. Guerra, 
M. Montenegro and one of the authors [10], as well as a simplification due to A. 
Moradifam [24]). It is also known that if N > 5, the problem has a unique singular 
radial solution for some A = As (as follows from the analysis in [2] and [10]) and, 
as in Barenblatt's result, that it has infinitely many regular radial solutions for the 
same value of the parameter (see the delicate work due to J. Davila, I. Flores and 
I. Guerra [11]). The case of general domains with Dirichlet boundary conditions 

■^For any solution of 111.211 . v is super harmonic. In particular, its spherical average v is a 
decreasing function of r. We will prove that v > 0, so that ?J(oo) > is always well-defined. 
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is essentially unexplored, due to the lack of a comparison principle. The case of 
homogeneous Navier boundary conditions, namely the equation 

j A 2 u = \e u in fi, 
I u = Au = on dfl, 



(1.8) 



where f2 is a smooth and bounded domain in R N , N > 1, seems, for now, easier to 
deal with. Our results continue a story initiated by C. Cowan, P. Esposito, and N. 
Ghoussoub in [TJ. They can be summarized as follows. 

Theorem 1.9. Let N > 1 and let Q be a smoothly bounded domain o/l w . Let u* 
be the extremal solution of U.S\ ). 

• If 1 < N < 12, then u* € C 00 ^). 

• If N > 13, i/ien u G C°°(r2 \ S) , where £ is a closed set whose Hausdorff 
dimension is bounded above by 

HdimiY.) < N - Ap* 
and p* > 3 is the largest root of the polynomial (X — i) 3 — 8(X — i) + 4. 

Remark 1.10. Theorem 11.91 was first proved for 1 < iV < 8 by C. Cowan, P. 
Esposito, and N. Ghoussoub, see [7]. As we were completing this work, we learnt 
that C. Cowan and N. Ghoussoub just improved their result to 1 < N < 10, see [8]. 
The question of partial regularity in large dimension was studied by K. Wang for 
the classical Gel'fand problem (see [31]). Unfortunately, we could not understand 
a part of his proof! an d our approach is somewhat different. We expect that the 
computed exponent p* is not optimal. Similarly, as observed by P. Esposito [33], 
the methods that we develop here will not yield the expected critical curve (resp. 
dimension) for the Lane-Emden system (resp. for the MEMS problem), for which 
new ideas must be found. It could be interesting to see if they improve the results 
obtained in the previously mentioned references [7J, [5], as well as the works of J. 
Wei and D. Ye [31], of J. Wei, X. Xu and W. Yang [3J], and of C. Cowan [B]. 



Finally, our Liouvillc-type result, Theorem [TTJJ will be used to prove the follow- 
ing. 

Theorem 1.11. Let 5 < N < 12 and let VI be a smoothly bounded domain ofW. N . 
Assume in addition that f2 is convex. Let u £ C 4 (f2) be any classical solution of 
li.<S|) and let v = — Au. There exists a compact subdomain lu a £1 such that if 

(1.9) / vdx< Kr N ~ 2 , 

JB r (x ) 

for every ball B r (xo) C co and for some constant K > 0, then, there exists a number 
C depending only on A, f2 ; N, K and the Morse index of u, such that 

IMU~ ( n) < C. 

Definition 1.12. In the above, the Morse index of a solution u is the number of 
negative eigenvalues of the linearized operator L = A 2 — \e u with domain D(L) = 
{u : u, Au £ H 2 (il)nHQ (il)}. According to standard spectral theory, this number 
is finite. 



the relation between the exponents in his interpolation inequality (3.18) seems incorrect. 
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Remark 1.13. If u is stable, then ()1.9j) automatically holds for some constant K 
depending only on fl,N, and u. See Lemma [5.21 We do not know whether this 
remains valid for solutions of bounded Morse index. In addition, it would be inter- 
esting to know how the number C depends on the Morse index of u. See the work 
of X.-F. Yang [37] for a result in this direction, in a subcritical setting. 

Notation. For any given function /, / connotes the spherical average of /. We 
write / < g (rcsp. / > g), when there exists a numerical constant C such that 
/ < C<? (rcsp. / > Cg). Br(x) denotes the ball centered at x and of radius R. For 
shorthand, Br = Br(0) and Ar is the annulus of radii R and 2R. 



2. Classification of stable radial solutions 

We prove here Proposition ll.il Take f3 > f3 and let u = up be the radial solution 
such that u(0) = and v(0) = j3. We claim that if /3 is large enough, then up is 
stable. We shall use the following inequality, found in 

(2.1) up<-^^r 2 forallr>0. 

Simply choose /3 > /?o such that 

?-3o 2 N 2 (N - A) 2 

(2.2) e-^ r < 1 , ' for all r > 0. 

16r 4 

Combining (|2.1[) . (|2.2[) , and the Hardy- Rellich inequality, we deduce that up is 
stable for /?>/?. So, we may define A = {(3 > {3q \ up is stable} and j3\ = inf A. 
By standard ODE theory, one easily proves that (3\ = min A. According to [3J, up 
is unstable. So, /3i > /3o- Also, by a result in p], solutions are ordered : if j3 > /?, 
then u^ < up. So, A is the interval [/3i, +oo). □ 

3. Construction of nonradial solutions 



We present here the proof of Theorem ll.3l Take a polynomial p of the form p.5p . 
Without loss if generality, we may assume that x° = 0. We look for a solution u of 
the form u = —p(x) + z(x), so that z and w = —Az satisfy 



-Az = w in R N , 

-Aw^erP^e* inR N . 



(3.1) 

For x e R N , let 

Z(x) = (1 + |^|2)2-iV/ 2 ( w{x) = {l + | a .|2 ) l-J V /2 > 

We claim that (^, IV) is a super-solution of (|3.ip . Indeed, straightforward calcula- 
tions yield 

-AZ > 2(N - A)W in R N , 

-AIV = Ar(Ar-2)(l + |x| 2 )- 1 ~ Ar / 2 in E^. 
Since Z < 1 and since we assumed that p(a;) > (1 + N/2)\x\ 2 in IR^, we have 

-AW > (1 + |^| 2 )- 1 - 7V/2 e^ > e-^/ 2 ^! V > e^^e 2 in R w , 
which proves our claim. 
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Since the system is cooperative, and (0,0) and (Z,W) form a pair of ordered 
sub- and super-solutions, we obtain the existence of a solution of (|3.1|) which further 
satisfies < z < Z and < w < W in M. N . Hence, u(x) := —p(x)+z(x) is a solution 
of O- 

To prove that u is stable outside a compact set, let us observe again that Z < 1 
in M^. So, we can find p > large such that 

(3.2) e u{x) < e z{x) e-^ < e 1 "^ < , M , for all Id > p. 

16|x| 4 

By the Hardy- Rellich inequality, u is stable outside _B p . Remark now that if 
min i= i...Af on > is large enough then (|3.2p is valid for all x £ R^, so u is stable in 
R N . This ends the proof of Theorem O □ 



4. Regularity of the extremal solution in dimension 1 < N < 12 

In this section, we prove the first part of Thcorcm ll.9l Let u denote the minimal 
solution to p.8[) associated to a parameter A G (A*/2, A*). Up to rescaling, we may 
assume that A = 1. The first ingredient in our proof is the following consequence 
of the stability inequality (jl.3l) : 



(4.1) / e%tp 2 dx< [ |Vv?| 2 dx for all <p e #o(fi). 

Jn Jn 

(|4.1j) was proved indepcndentljfl by C. Cowan and N. Ghoussoub in [5] and A. 
Farina, B. Sirakov, and one of the authors in |13j . See e.g. Lemma 16.11 below for 
its proof. Now, Let us write the problem (|1.8p as a system in the following way: 

{—Au = v in 17, 

-Av = e u in O, 

u = v = on dfl. 

Fix a > | and multiply the first equation in (|4.2[) by e au — 1. Integrating over fi, 
we obtain 

/ (e au -l)v dx = a [ e au \Vu\ 2 dx = - f | V (e^ - l) | 2 dx. 
Jn Jn a Jn 



By 

/ eT ( e ^r _ i) 2 dx < \W (e^r - l)| 2 dx. 
in in 
Combining these two inequalities, we deduce that 

(4.3) f e% (e^ - l) 2 dx < 7 / (e QM - l)v dx. 

Similarly, multiply the second equation in (|4.2j) by i> 2q_1 and use (|4.ip to deduce 
that 

(4.4) / e*v 2Q dx < — / e u v 2a ~ 1 dx. 

V ' Jn ~ 2a - 1 In 



^The result was first made available online by C. Cowan and N. Ghoussoub. See also the 
Appendix for another result found by several authors. 



THE GEL'FAND PROBLEM FOR THE BIHARMONIC OPERATOR 



7 



By Holder's inequality, 



eV"- 1 dx<( e%v 2a dx) / e*e QU dx ) and 

2a-l 

*u 2a dx) II e*e Q " dx 
n 



Plugging these inequalities in (|4.4p and (|4.3[) respectively, we deduce that 

.,2 



a 

2^ 1 



*v 2a dx) <- -[ I e^e au dx) and 



e f - l) 2 dx<jyj e%v 2a dx 



Multiplying these inequalities, it follows that 

f e * (e* - l) 2 dx < -2— [ e {a+ ^ u dx 

so that 

( 1 j / e%e au dx<2 [ e^ u dx. 

Apply Holder's inequality to the right-hand side. Then, 

Q + l 

and so 

(4.5) (l--£J)(J a .-P. 

Let a* > 5/2 denote the largest root of the polynomial X 3 — 8X + 4. We have 
just proved that e" is uniformly bounded in L p (il) for every p < p* = a* + |. In 
particular, e u is bounded in L P (Q) for some p > 3. Using elliptic regularity applied 
to (|1.8j) . this implies that u* is bounded, hence smooth whenever TV < 12. □ 



5. Growth of the L 1 norm 



This section and the next provide preparatory results that will be used both for 
the proof of the Liouville-type theorem and the partial regularity result. We begin 
with the case of R N . 

Lemma 5.1. Assume N > 5 and let u be a solution of (|1 .2|) which is stable (resp. 
stable outside the ball Br ). Let v = —Am, v its spherical average, and assume 
that v(oo) = 0. Let Br denote the ball of radius R (resp. the annulus of radii R 
and 2R). Then, there exists a constant C depending only on N (resp. on N,u,Rq) 
such that 

(5.1) / v dx < CR N ~ 2 for every R>0 (resp. R > R„). 
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Proof. We claim that v > 0. Since the equation is invariant under translation, it 
suffices to prove that v(0) > 0. Let v be the spherical average of v and assume by 
contradiction that u(0) = v(0) < 0. We have 



(5.2) 



-Av = e u in t^, 
-Au = v mR N . 



In particular, —Av = — r 1_w (r A ' _1 v')' > 0, and so v is a decreasing function of 
r > 0. Since we assumed that v(0) < 0, it follows that v{r) < for all r > 0. So, u 
is subharmonic. It follows that u is an increasing function of r > 0. In particular, 
u is bounded below and so, given R > 2R$, 

(5.3) f e"dx>e"(°) f dx > R N . 

J B2r\Br J B2r\Bh 

Take a cut-off function ip e C£ (B 4 \ Bi), with < ip < 1 in R w and tp = 1 in 
B2\B\. Applying Jensen's inequality and the stability inequality (|1.3p with test 
function <p(x/R), we find 

e"dx< j ~^dx<R N - A 

which contradicts (|5.3j) . if i? is chosen large enough. Hence, v(0) > 0. Now, if u is 
stable, take r > 0, take a standard cut-off function ip £ C^iB-i) such that < (p < 1 
and (/3 = 1 in £?i. Apply the stability inequality (|1.3j) with test function ip(x/r) to 
get 

JV-4 



(5.4) / e u < 



If u is only stable outside Br , fix r > 2i?o and take a cut-off function <p <G (R ) 
such that ip = in -Br , <p = 1 in \ Br _|_i, <p = outside -B2r, and |A</?| < r~ 2 
in £?2r \ ■ Using the stability inequality (|1.3[) with this test function we find that 



(5.5) / e" dx = / e u dx + / e" dx < 1 

J B r J Br q + i J B r \Bji Q .\-i 

Now, we rewrite the first equation in ([5.2ft as 

Integrate on (0,r). By (|5.4p . which holds for every r > if u is stable (resp. by 
()5.5|) . which holds for r > 2R,q if u is stable outside Bn ), 

r r 



-r n - L v'(r) = / t^-^dt^r 1 
Jo 

We integrate once more between R and +00. Since TJ(oo) = 0, we obtain 

v(R) < R- 2 . 

Clearly, dET} follows. □ 



Now, we turn to an analogous result on bounded domains. Let us recall that the 
extremal solution is stable. 
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Lemma 5.2. Let N > 3 and let u be the extremal solution of U.8\) . Set v = — Au. 

Fix io 6 ^ and let Rq = dist(xo, dCl)/2. Then, there exists a constant C > 
depending only on N , £1, and Rq such that for all r < Rq, 

(5.6) f vdx< Cr N -' 2 . 

JB r {x a ) 

Proof. Without loss of generality, we may assume that xq = and A* = 1 . 
Step 1. There exists a constant C = C(N, Rq) such that for all r € (0, Rq), 

(5.7) f e u \x\ 2 ~ N dx < Cr~ 2 . 

J[r<\x\<R ] 

To see this, consider the function ip '■ ^ir q — > R N given by 



(5.8) iP(x) 



a - b\x\ 2 if |x| < r, 



| x |l-W/2 if r < | x | < 2i? , 

where the constants a = K±l r l - N / 2 anc i 5 = N=2 r -i-N/2 are choggn so i s 
C 1 and piecewise C 2 . Take a standard cut-off function 

(5.9) (eC 2 (B 2 ) such that < C < 1 and C = 1 on Bi, 

so that (p(x) = ip(x)C(x/Ro) belongs to H&(Q) n H 2 (n). 
Since u is stable, we have 

(5.10) f e"|aj| a - w dx < [ e u tp 2 dx < [ \A<p\ 2 dx. 
J[r<\x\<R ] Jn Jn 

By 

(5.11) f \Aip\ 2 dx < 4,N 2 b 2 [ dx< Cb 2 r N < C'r- 2 . 

J[\x\<r] J[\x\<r] 

Similarly, 

(5.12) f \Acp\ 2 dx<C [ 

J[r<\x\<R ] J[r 

Finally, 

(5.13) f \Aip\ 2 dx < C(N, R ) < C'r' 2 . 
J[R a <\x\]nn 

Collecting ([5U0| . (|5TT|l . (f57T2jl . and (|5~T3j) . the estimate (J5J]) follows. 
Step 2. Take as before a standard cut-off £ satisfying (|5.9[) and let 



\x\~ 2 ~ N dx < C'r 

r<|s|] 



V(x)=/ \x-y\ 2 - N ((y)dy, for x € R N . 

Then, there exists a constant C, depending on N only such that 
(5.14) ip(x) < C*min(l, \x\ 2 ~ N ), for 1 e M". 

Indeed, if |x| > 4 and \y\ < 2, 

\x-y\ > \x\ - \y\ > ^\x\ 

so that 

^x)<2 N ' 2 \x\ 2 ~ N f C(V) dy < C\x\ 2 ~ N , 
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while if |x| < 4 and \y\ < 2, 

\x-y\ < \x\ + \y\ < 6 

and so 

lK*0 < HCIk-(R*) / \x - y\ 2 - N dy < C. 

J\x-y\<6 

(f5T4|) follows. 

Step 3. Take r £ (0, Rq), £ a standard cut-off function satisfying (|5.9[) . and let ip r 
be the solution to 

| — Aif r = Q(x/r) in £7, 
) (p r = on 9f2. 



(5.15) 



Then, there exists a constant C depending on ./V only such that 

(5.16) ip r (x) <Cr 2 mm (l,r N - 2 \x\ 2 - N ) , for all x e ft. 

This easily follows from the maximum principle, observing that a constant multiple 
of r 2 ip(x/r) is a supersolution to the above equation. 

Step 4. There exists a constant C depending on N and fl only, such that 



(5.17) / e u dx<C. 

This is an obvious consequence of Equation (|4.5[) (which holds for any a £ (^,5/2]) 
and Holder's inequality. 

Step 5. Multiply (|5.15[) by v = —Au and integrate by parts. Then, 

v dx < / v(,(x/r)dx= I e u ip r dx. 

x\<r] Jn J Q 



Using Step 3, we have 



f e u ip r dx < Cr 2 [ e u dx. 

J\\x\<r] J\\x\<r] 



[|x|<r] J[\x\<r] 

Using stability with test function £(x/r), we also have 

e u dx < Cr N ~ A 

x\<r] 

and so 

„N-2 



f e V dx < Cr 1 



By Step 1 and Step 3, 



f e u ip r dx < Cr N f e u \x\ 2 - N dx < Cr N - 2 . 

J[r<\x\<Ra\ J[r<\x\<R a ] 

Finally, using ([BTTT) . 

f e u <p r dx < CR 2 Q - N r N f e u dx < C'r N ~ 2 . 

(EH) follows. □ 
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6. A BOOTSTRAP ARGUMENT 

Our next task consists in improving the L^-estimates of the previous section to 
L p -estimates for larger values of p. We do this through a bootstrap argument which 
is reminescent of the classical Moser iteration method, up to one major difference: 
we will take advantage of both the standard Sobolev inequality and the stability 
inequality (|1.3[) . To be more precise, rather than (|1.3p . the following interpolated 
version of it will be used. 

Lemma 6.1. Let fl be an open set of M. N , N > 1. Assume that the stability 
inequality (| 1 .3|) holds. Then, for every s g [0, 1], 

/ e su <p 2 dx< [ \(-A) s / 2 <p\ 2 dx for all ip e H 2 (Q) (1 H^(Q). 
Jn Jn 

In particular, for s = |, 

(6.1) / e^ip 2 dx < I \\7ip\ 2 dx for all ip g Hq(Q). 

Jn Jn 

Proof Consider first the case where O = M. . We apply complex interpolation 
between the family of spaces X s , Y s given for < s < 1 by 

X s = L 2 ((2nr N \^ s dO , Y s = L 2 (e su dx). 

Recall that the inverse Fourier transform J 7 " 1 : X a —> Y satisfies ||^||jC(x ,Yb) = ^" 
Furthermore, by the stability inequality (|1.3|) and Plancherel's theorem, we have 

/ e u ip 2 dx< f \Aip\ 2 dx = (2tt)-" f \£\ 4 \\T((f>)\ 2 d£ 
Jr n Jm n Jwl n 

Thus, J 7 ' 1 : Xi ->• Yx satisfies || TWax^Y,) < 1- By [H Theorem 2], we deduce 
that W^W C (x s ,Y s ) < 1 for all < s < 1. 

In the case where f2 is a bounded open set, simply repeat the above proof, using 
the spectral decomposition of the Laplace operator in place of the Fourier transform. 
More precisely, let (A&) denote the eigenvalues of the Laplace operator (with domain 
D = H 2 (fl) n Hq(Q)), let (pk be the fc-th component of <p in the Hilbert basis of 
cigcnfunctions associated to (Afe) and interpolate between the family of weighted 
L 2 -spaces X s , Y s corresponding to the norms 

IM&.=£ A fc*l&| a ' 11*.= [ e™<p 2 dx. 

In the case where f2 is an unbounded proper open set, take k > 0, let = £1 n B k 
and 

- Mfc = inf|^ (|A<^| 2 - e > 2 ) dx : v eH 2 (n k )nHi(n k ),M LHnk) = iy 

Then, the previous analysis leads to 

/ [{e u - Vk)+Y P 2 dx < f \{-A) s / 2 ip\ 2 dx for ail <p e H 2 {tt k ) n H%(n k ). 
Jn k Jn k 

By (|1.3| . lim/o^+oo — fik > and the result follows. □ 
Our next lemma is simply the first step in the Moser iteration method: we 
multiply the equation by a power of its right-hand side, localize, and integrate. 
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Lemma 6.2. Let fl be an open set ofR N , N > 1. Assume (u,v) <E C 2 (f2) 2 solves 

-Au = v in f2, 

-Av = e u in £1. 

Take a > |, ip £ (7^(0). Then, there exists a constant C depending on a only, 
such that 



(6.2) V2a 1 \\V(v a <p)\\ L 2 m < \\e%v a -^\\ L 2 {n) + C\\v a \7 v\\mn)- 



and 
(6.3) 



|V(e*»|| La(n ) < \\e% u v*tp\\ L 2 iu) +C\\e% u V<p\\ L 2 {n) . 



Proof. Since the computations are very similar, we prove only (|6.2[) . Wc multiply 
—Av = e v by v 2a ~ 1 Lp 1 and we integrate. We obtain 



e u v 2a ~ V dx = / Vv ■ V(v Za - L ip 2 ) dx 
Jn 

2a - 1 , ,_,„ |2 ,„2 



or 
2a - 1 



Of 

2{a-\) 



W*| V dx + 2 v Za - L ipVv ■ Vip dx 



|V(v»| 2 dx - / v 2a \Wif\ 2 dx 



v a <jN<p ■ Vv a dx. 



In the last term of the right hand side, replace tp\7v a by V(v a ip) — v a V<p. Then, 



eV Q " V dx = 



2a -I 

a 2 Jn 
2(a-l) 



V(«»r dx- 

a Jn 



v a V V V(v a ip) dx, 



v 2a \Vip\ 2 dx 



which we rewrite as 

{2a -1) / |V(u»| 2 dx = 2(a - 1) / v a \7 tp\7 (v a tp) dx+ 



(6.4) 



By the Cauchy-Schwarz inequality, 



v 2a \Vip\ 2 dx + a 2 / e u v 2a ~ L ip 2 dx 



v a V<pV{v a tp) dx 



< 



|V(v»| 2 dxj (j v 2a \V<p\ 2 dx 

Plugging this in (|6.4[) , we obtain a quadratic inequality of the form 
{2a - l)X 2 < 2\a - 1\AX + A 2 + B 2 

where 

X = ||V(t>»|| 2 , A=||« a Vvj|| a and B = a\\eMv a -^ V \\ 2 
Solving the quadratic inequality, we deduce that 



x < W - M + Vh jg + (2a 1){A 2 + B 2 ) < B +C A. 



2a- 1 



V2a- 1 



THE GEL'FAND PROBLEM FOR THE BIHARMONIC OPERATOR 



13 



(|6.2|) follows by replacing A, B and X with their values. □ 
We have just used the equation. Now, we use the stability assumption. 



Lemma 6.3. Make the same assumptions as in Lemma w.SX Assume in addition 
that (|6.1[) holds for every f £ C*(f2). Let a", a* denote the largest two roots of the 
polynomial X 3 — 8X + 4. Then, for every a £ (a", a*), there exists a constant C 
depending on a only such that 



(6.5) / \V(v a tp)\ 2 dx < C I v 2a \V<p\ 2 dx 

n 



or 



(6.6) f |V(et»| 2 dx < C [ e au \\7tp\ 2 dx, 

Jn Jq 

for any ip £ Cjr(f2). 

Proof. By Holder's inequality, 



e u v Za - L tp z dx < [ e" V dx) [ e^v 2a ip 2 dx 
Using the stability inequality (|6.1j) . we deduce that 

eV Q -V dx<HiK 2 ~i, 



where we set 
Similarly, 



F=||V(e*»|| 2 and K = ||V(w»|| 2 . 



e au vtp 2 dx < K»H 2 ~ 

i 

Combining with (|6.2p - (|6.3p . this gives 



(6.7) V2a 1 H < K^H 1 -^ +C||e^"Vy>|| 2 , 

a 

(6.8) ^=K < H^K 1 -^ +C\\v a V<p\\ 2 . 
Multiply dUTJ) by ([5T5]) . Then, 



(6.9) ( 2 ^ 2a 1 _ i ) HK < aH^K 1 -^ + bK^H 1 "^ + ab, 



where 

a = C\\e% u X7tp\\ 2 and b = C\\v a \7tp\\ 2 . 
Note that for a £ (a", a*), 



Introduce 



2V2a - 1 
8 := — — = 1 > 0. 



X = K^H 1 -^ and Y = H^K 1 -^. 



Then, (|6 .9[) can be rewritten as 

5XY <aY + bX + ab, 
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and so, either X is bounded by a multiple of o or Y by a multiple of b. In the former 
case, recalling (|6.7[) . we obtain (|6.6|) . In the latter case, ()6.8[) implies ()6.5|) . □ 



In the two previous lemmata, we have used successively the equation and the 
stability assumption. Now, we apply the Sobolev inequality to set up a bootstrap 
procedure. 



N-4a 



Lemma 6.4. Make the same assumptions as in Lemma \ 6.3[ Take a £ (a%a*) 
and for R > 0, let Br denote a ball of radius R (resp. an annulus of radii R and 
R/2) contained in fl. Assume that there exists a constant C depending on N and 
a only, such that for all R > (resp. for all R large enough) 

(H a ) [ (e au + v 2a ) dx<CR 

J B R 

Then, ( H n „. J also holds. 

\ N-2 UL J 

Bootstrapping the above lemma, we find 

Corollary 6.5. Make the same assumptions as in Lemma \6.3[ Assume that §H, 
holds for some a <G (a", a*). Then, 



[ e pu dx < CR N - ip , for all p < p* := a* + -, 

J Br 2 

r ON 

/ v q dx < CR N ~ 2q , for all q < q* := a* 

J Br ~ N-2 



(6.10) 



Remark 6.6. The inequality §H a l can be further simplified if the boundary val- 
ues/limiting behavior at infinity of the solution is known. See in particular Propo- 
sition IA.1I 



Proof of the lemma \6^\ Assume \H a \ is valid. By Lemma [B~3l either (|6.5I) or (|6.6p 
holds. 

Assume that (|6.5|) is valid (the other case is similar). Using the Sobolev embedding, 
we obtain 



Ta 2* 



ip 2 dx) < \V(v a ip)\ 2 dx< / v 2a \V<p\ 2 dx 



Take a standard cut-off function ip <G C^B^) such that O<-0<1, -0 = 1 in B\, 
and = outside B 2 . Apply the above inequality with tp(x) = ip(x/R) and use 
TO- Then, 



(6.11) / v Ta dx < R N - 2 Ta . 

J Br 

Going back to ()6.7|) . we deduce similarly that 

(6.12) f e^ u dx<R N - 2 '- 2a . 



Br 

□ 



Proof of the corollary 1 6. 51 By Holder's inequality, if Qff a | ) holds for some a, then 
{Hp) holds for all (3 < a. So, bootstrapping jH a l , we easily deduce that it holds 



for all a < jj^_ 2 a* . Fix at last a < a* and apply now stability (|6.1[) . with test 
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function e^ u ilj(x/R), to deduce that (|6.10[) holds for all p = a + | < p* = a* + h. 

□ 



7. The Liouville theorem 

We prove here Theorem ll.7l Assume by contradiction that there exists a solution 
u of (|1.2[) which is stable outside a compact set K C Br and such that TJ(oo) = 0. 
Step 1. e u e L P (R N ), for every p < p* . 

By Lemma 15. 11 we have 



(7.1) / v dx < CR N - 2 for every R > R , 

A,. 



where An is the annulus of radii R and 2R. In addition, stability (|1.3[) implies 
that 



(7.2) 



/ e u dx < CR N ~ 4 for every R > R a 



Recall now the following standard elliptic estimate : for p G [1, jfz^), 

IM|z,f>(B 2 \-Bi) ^ l|Af|| L i( B2 \ Bl ) + |MUi(B 4 \B 1/2 ), 

and its rcscalcd version 

R~%\\v\\l*(a r ) < R~ N (i? 2 ||Av|| L i (Afi) + \\v\\ L i iBiR \B R/2 ) 

Applying this estimate respectively to u and v, we deduce from (|7.1I) . (|7.2j) that 
( |H a D holds for any a <S [1, 7735) and all large R. Hence, (|6.10[) holds. By a 
straightforward covering argument, Step 1 follows. 

The rest of the proof is very similar to the one given in [9]; 
Step 2. 

lim \x\ 4 e u = 0. 

\x\— >+oo 

By Step 1., given any S > 0, we can choose R large enough such that 



(7.3) / e pu dx < 5 P . 

J\y\>R 

Let x G R n \~Br^, \x\ > 4R. Set R = This yields 

B{x, R/4) C A fl = B 2R \B R c {y G R N I |y| > i?}. 

Thus, we have 

/ e pu dx <CR N - 4p , for all p < p* := a* + \, 

f74l ^B(x,fl/4) 2 

/ da; < CR N ~ 2q , for all a' 

JB(x,R/4) N - 2 

Next fix (5 > and consider w = e u . By Kato's inequality, w satisfies 
(7.5) — Aw — vw < in R N . 
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Take e small enough such that t^tj < wz2 a * ■ Here, we have used the assumption 
5 < N < 12. The Serrin-Trudinger inequality [571[3D| for subsolutions to ([73)1 
ensures that for any p < p* 

IMlL-(B(z,f )) < CR P NliP(B(x,# )) 

where C depends to N, p and 

In particular, for p = ^ and using (|7.4I) 
(7.6) e"^<Cir 4 |H| L „ /4(s(:Cif)) 
where C depends only to AT. Combining ()7.3[) and (|7.6|) gives 

e u(«) < C5R- i < Six]' 4 , 

which proves Step 2. 

Step 3. By Step 2., there exists Ri > Rq such that 

-Av < —7 for all r > R\. 

~ 2r 4 

Hence, 

^) > ~ 2(jV 1 4)r3 forallr>i?, 

Integrating between r and +oo, this yields 

2^2 ^2(AT-4) ~ 7 
Since N > 5 and choosing i?2 > -Ri large enough to have 

1 C'(N) 



Au(r) = v(r) = v(r) — v(oo) < ^-^ ( 777^7 — „jy_ 4 ) ^ or a ^ r > 



2(iV-4) r w - 4 

we get 



< 1 for all r > R 3 , 



-Au < — ^7 for all r > R 3 . 
2r 2 



In the same way, we have for some R 3 > R2 

u'(r) >— for all r > R 3 . 

r 

Integrating the latter and taking the exponential, we get 

r 4 e «(r) > cr 3 for all r>R3i 

where the constant c > docs not depend on r. By Jensen's inequality, we have 
for any r > R 3 , 



or 



3 < r i e u(r) < r 4g?T^ < max |a;|4 e f(*). 



x — r 



This contradicts Step 2. □ 
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8. Partial regularity of the extremal solution in dimension TV > 13. 

In this section, we prove the second part of Theorem 11.91 As in the previous 
sections, we interpret our equation as the system (|4.2p . Here, our task consists in 
showing that if a rescaled LP norm of e" is small on some ball B, then it remains 
small on any ball of smaller radius, which is included in B. This provides an 
estimate in Morrey spaces, for which an e-regularity theorem is available, thanks 
to the Moser-Trudinger inequality. 

Let u be the extremal solution of (|1.8[) . By scaling, we may assume that 
■Bi(O) CC ft. For any x £ Bi (0), 1 < p < p* and < r < 1 - \x\ we define 

E(x,r)=r 4p ~ N f e pu dy. 

J B r (x) 

Lemma 8.1. Assume E(0, 1) < 1 and fix p < q < p* . Then, there exists a positive 
constant C = C(7V,p, q) > such that 

i p 1 

(8.1) E{x,r) < CE(0, I) 1- " for all x G Bi(0), < r < -. 

2 2 

Proof. Let u = u\ + 112, where 

j-A Ul =0 inBi(O), j-Au 2 =V in Bi(0), 

j ui=u on dB 1 (0), j u 2 = on 9Bi(0). 

By the maximum principle, Ui,U2 > in 73i(0). Write i?(x, r) = £1 + i?2, where 

El = r i P -N f p pu dy ) ^ = r 4 P -Ar /" e p« 

JB r (.T)n[u 2 <e] JB r (i)n[ii 2 >8] 

and where 9 > will be fixed later on. By Kato's inequality, e pui is subharmonic. 
By the mean value inequality, it follows that 

(8.2) i e pui dy<£ e pui dy. 

J B r {x) JBi{x) 
2 

Hence, 

1.3) E l < r 4p - N e pf> [ e pui dy < r Ap e p!> [ e pu dy 

JB r (x) JB!(x) 

< r ip e pf> [ e pu dy < r 4p e p9 E(0, 1). 

7Bi(0) 



To estimate E2, we first note that 

(8.4) \B r (x) R [u 2 >6}\< e- pS f e pu dy < ^f^' - 

J B r (x)n[u 2 >e] r ip e p » 

By Holder's inequality, we also have 

(8.5) / e pu dy<([ e qu dy) {B r (x) n [u 2 > 6}} 

J B r {x)n[u 2 >9] \JB r (x) J 

Recall from Section [6] that 

e gu dy < r N~± q ^ 

B r {x) 
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Using this together with (|8~5j) . we find 



.6) E 2 < r ip ~ N 



e qu dy 



B r (x) 



{B r (x)n[u 2 > 6]} 



< 



E{x,r) 



1-E 



Combining (|E31) and ([EI]), for all x e Bi (0) and < r < \ we find 



A/', 



(8.7) E(x, r) < r ip e 9p E(0, I ) 4 

If for some x £ Bi (0) and < r < i we have 

~E(x, r y ] 



E{x,r) 



2 

1-E 



- 4p £(0, 1) > 



,.4/) 



then, using the assumption E(0, 1) < 1, we find 

2g — p 



E(x,r) < r 4p ^E{0,l)^ < ^(O,!) 1 "?, 



which is the desired inequality. So, assume that (|8.8j) docs not hold. Then, we can 
find 6 > such that 

i- £ 

r**{<?°f-*E(Q,l) 



E(x,r) 



E(x,r) 



r 4p 

that is, 

r 4p e p6 E{0,l) 
Using this fact in (|8.7j) we derive 

E(x, r) < r 4p e 8p E(0, 1) < \e(0, 1)E(x, r) 
and finally E(x,r) < E(0, 1) 1 "! 



J7 P_ 

2<j-p 



□ 



Proof of Theorem 11.91 completed. Let £ be the singular set of u. Let 1 < p < 

q < p* be fixed. We claim that 



(8.9) 



Sc uefi: lim sup r 

r->0 



4p-N 



e pu > 



B r {x) 



Assume to the contrary that there exists xo £ £ such that lim r _).o E(xq, r) = 0. Fix 
£o > 0. Then, there exists p = p(eo) > such that E(xo 7 r) < so for all < r < p. 
For the sake of clarity let us assume that xq = and p = 1 . Thus, we can apply 
Lemma 18. f I and obtain 

l-E 1 

E(x,r) < C(N,p,q)e q for all x e Bi (0) , < r < — . 

2 2 

Using Holder's inequality this implies that e u € M N / 4 (Bi (0)) and 
\\e u \\ MN/HBim <C(N,p,q)el 



1_ 1 

. V <3 



Let 



w(x) 



\x-y 



4-N u{y) 



dy, xeBi{0) 



THE GEL'FAND PROBLEM FOR THE BIHARMONIC OPERATOR 



1!) 



and w = u - w. By Lemma 7.20 in [20] it follows that e Pw € L x {Bi{Q)) for all 

C{N) 



< p < 



e" jv 

"mt fli o 



Since w is biharmonic in Bi (0), it follows that e" € L^ 3 (.81/4(0)) for all 

0<^<C(JV,p,g)4~*- 
Letting eo << 1 small, we have e u € L°(B 1 u(0)) for some /3 > N/4, which, by 
standard regularity theory, yields u G L°° (B 1 / 8 (0)) and contradicts 6 E. 

Hence, (|8.9p holds for all 1 < p < p* . By [HI Lemma 5.3.4], it follows that 
■H w - 4p (£) = for all 1< p < p* . Thus U dim (£) < N — Ap* . □ 

9. A SCALING ARGUMENT 

This last section is devoted to the proof of Theorem II .111 By rescaling, we may 
always assume that A = 1. By standard elliptic regularity, it suffices to show that 
u < C in D,. We assume to the contrary that there exists a sequence of solutions 
u n of fixed index k, such that M n := max^u„ — > +oo, asm +oo. Let x n denote 
a corresponding point of maximum of u n . Passing to a subsequence if necessary, 
we may assume that there exists xq € fl, such that x n — > xo, as n — >• +oo. Since f2 
is convex, we can also assert that u„ is uniformly bounded in a fixed neighborhood 
of the boundary i.e. x n € ui CC f2 for large n. See e.g. (33] for this standard 
boundary estimate. 

We use a scaling argument. Let r„ = e~ M "/ 4 and U n (x) = u n (x n + r n x) — M n , 
for a; e fi„ := i(Q — x„). Then, [/„ solves 



(9.1) 



A 2 L7„ = e Un in O 



U n + M n = AU n = on dVL n . 

Since x n <G u> CC 0, fi„ — > R , as n — > +oo. We claim that (U n ) is uniformly 
bounded on compact sets of R . To see this, fix a ball Br and n so large that 
B R c n n . Write -AU n = V n = V n l - V%, where V* solves 

f -AV; 1 = e 17 " in B R , 

\ V,l = on dB R , 

and is harmonic in Br. Since U n < 0, V 4 is uniformly bounded in Br. Using 
the assumption (|1.9[) , we also have 

N-2 



(9.2) / «„ dx < C{r n R) 

JB{x n ,r n R) 

In other words, V n = —AU n is bounded in ^(Br). Since is uniformly bounded 
in Br, it follows that is bounded in ^(Br). Since V% is harmonic, is 
uniformly bounded in Br/ 2 . Similarly, write U n = U„ — U%, where 

f - AC/ 4 = V n in B R/2 , 

\ 1/1=0 on dB R/2 , 

and U% > is harmonic in Brj 2 . Since V n is bounded in Br/ 2 , so is ?7„. Since 
Un{0) = U^(0), we may now apply Harnack's inequality to conclude that U% is 
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uniformly bounded in -Br/4. Hence, U n is uniformly bounded on compact subsets 
of M. N . So, we may pass to the limit in the first line of (|9.1[) and find a solution U 
of finite Morse index to ([L~2)l. Thanks to flgHJ}, if V = -AU and R > 0, then 

f V dx < CR N ~ 2 , 
JB R 

which is possible only if V{oo) = 0. This contradicts Theorem 1 1.71 

□ 



Appendix A. Appendix 

The following is a consequence of Kato's inequality. It was also proved by N. 
Ghoussoub et al. [19] in unpublished work. Its proof builds upon a similar result 
of P. Souplet |28j . as well as similar inequalities in bounded domains obtained by 
C. Cowan, P. Esposito, and N. Ghoussoub [7]. 

Proposition A.l. Let u be a stable solution of \l. e 3!) and v = — Au. Then, 

v > V2e% in R N . 

Proof. Let w := \[2e^ — v. A straightforward calculation yields 

(A.l) Aw > -^-=e3w in R N . 

Since v > (see the proof of Lemma l57Tj) . we have w + < \f2e% in R N . We multiply 
(|A.1|) by w + and integrate over Br(0), for R > 0. We obtain 

(A.2) f |V U >+| 2 + ^ f e^\w+\ 2 < [ w + — 



V2 J J dR 

Br(0) B r (0) dB R (0) 



For all R > define 



Then (fA3|) reads 



f (£)'•= I J \w + (Rx)\ 2 da(x) 

dBi(0) 



(A.3) j \Vw + \ 2 + ^= J e^\w+\ 2 <R N - 1 f'(R). 

B R (0) Br(0) 

Since u is stable and e u > 2\w+\ 2 in R N , we have 

R N ~ 4 > J e u ds> J r N - l f{r)dr. 

Br(0) 

In particular, / cannot be strictly increasing in a given interval [S, +oo). Hence, 
there exists an increasing sequence {Rj} such that Rj — > oo and f'(Rj) < 0. Now, 
letting R = Rj in ([A3]) , we find w + = 0, that is, v > \f2e% in R^. 

□ 
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